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ABSTRACT 

Apart from omnidirectional, a solid elastic sphere is a natural multimode and mul- 
tifrequency device for the detection of Gravitational Waves (GW). Motion sensing 
in a spherical GW detector thus requires a multiple set of transducers attached to 
it at suitable locations. Resonant transducers exert a significant back action on the 
larger sphere, and as a consequence the joint dynamics of the entire system has to 
be properly understood before reliable conclusions can be drawn from its readout. In 
this paper, we present and develop a mathematical formalism to analyse such dynam- 
ics, which generalises and enhances currently existing ones, and which clarifies their 
actual range of validity, thereby shedding light into the physics of the detector. In 
addition, the new formalism has enabled us to discover a new resonator layout (we 
call it PHC) which only requires five resonators per quadrupole mode sensed, and 
has mode channels, i.e., linear combinations of the transducers' readouts which are 
directly proportional to the GW amplitudes. The perturbative nature of our proposed 
approach makes it also very well adapted to systematically assess the consequences of 
small mistunings in the device parameters by robust analytic methods. These prove 
to be very powerful, not only theoretically but also practically: confrontation of our 
model's predictions with the experimental data of the LSTJ prototype detector TIGA 
reveals an agreement between both consistently reaching the fourth decimal place. 
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1 INTRODUCTION 

The idea of using a solid elastic sphere as a gravitati onal wave (GW ) antenna is almost as old as that of using cylindrical 
bars: as far back as 1971 Forward published a paper ( Forward 1971 ) in which he assessed some of the potentialities offered 
by a spherical solid for that purpose. It was however Weber's ongoing philosophy and practice of using bars which eventually 
prevailed and developped up to the present date, with the highly sophisticated and sensitive ultracryogenic systems currently 
in operation — see (Coccia et al. 1995b) and (Coccia et al. 1997) for reviews and bibliography. With few exceptions (Ashby & 



Dreitlein 1975; Wagoner & Paik 1977), spherical detectors fell into oblivion for years, but interest in them strongly re-emerged 



in the early 1990's, and an important number of research articles have been published since which address a wide variety of 
problems in GW spherical detector science. At the same time, international collaboration has intensified, and prospects for 
the actual construction of large spherical GW observatories (in the range of ~100 tons) are being currently considered in 
several countries Q even in a variant hollow shape (Coccia et al. 1998). 

A spherical antenna is obviously omnidirectional but, most important, it is also a natural multimode device, i.e., when 
suitably monitored, it can generate information on all the GW amplitudes and incidence direction (Magalhaes et al. 1995), 



a capability which possesses no other individual GW detector, whether resonant or interferometric (Dhurandhar & Tinto 



198£). Furthermore, a spherical antenna could also reveal the eventual existence of monopole gravitational radiation, or set 



thresholds on it (Bianchi et al. 1998). The theoretical explanation of these facts is to be found in the unique matching between 



* There are collaborations in Brazil, Holland, Italy and Spain. 
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the GW amplitude structure and that of the sphere oscillation eigenmodes: a general metric GW generates a Udal field of 
forces in an elastic body which is given in terms of the "electric" components Roioj (t) of the Riemann tensor at its centre of 



mass by the following formula (Lobo 1995) 



fGw(x,t) = 



(1) 



1=0 and 2 
m— — L, . .. ,1 



where f'''"'^(x) are "tidal form factors", while are specific linear combinations of the Riemann tensor components 

RoiOjit) which carry all the dynamical information on the GW's monopole {l~0) and quadrupole (1 = 2) amplitudes. It is 
precisely these amplitudes which a GW detector is aimed to measure. 

On the other hand, a free elastic sphere has two families of oscillation eigenmodes, so called toroidal and spheroidal 
modes, and modes wit hin either f amily group into ascending series of /-pole harmonics, each of whose frequencies is (2/+1)- 
fold degenerate — see ( Lobo 1995 ) for full detail s. It so happens th at only monopole and/or quadrupole spheroidal modes can 
possibly be excited by an incoming metric GW ( Bianchi et al. 1996 ), and their GW driven amplitudes are directly proportional 
to the wave amplitud es ff^'"'''(t) o f equation (0). It is this very fact which makes of the spherical detector such a natural one 



for GW observations ( Lobo 1995 ). In addition, a spherical antenna has a significantly higher absorption cross s ection than a 
cylinder of like fundamental frequency, and also presents good sensitivity at the second quadrupole harmonic ( Coccia et al. 
1995a|. 

In order to monitor the GW induced deformations of the sphere motion sensors are required. In cylindrical bars, current 
state of the art technology is based upon resonant transducers (Astone et al. 1993; Hamilton et al. 198£| ) . A resonant transducer 
consists in a small (compared to the bar) mechanical device possessing a resonance frequency accurately tuned to that of the 
cylinder. This frequency matching causes back-and- forth resonant energy transfer between the two bodies (bar and resonator) , 
which results in turn in mechanically amplified oscillations of the smaller resonator. The philosophy of using resonators for 
motion sensing is directly transplantable to a spherical detector — only a multiple set rather than a single resonator is required 
if its potential capabilities as a multimode system are to be exploited to satisfaction. 

A most immediate question in a multiple motion sensor system is: where should the sensors be? The answer to this 
basic question naturally depends on design and purpose criteria. Merkowitz and Johnson (M&J) made a very appealing 
proposal consisting in a set of 6 identical resonators coupling to the radial motions of the sph ere's surface, and occupying 
the positions of the centres of the 6 non-parallel pentagonal faces of a truncated icosahedron ( Johnson fc Merkowitz 199^ ; 



Merkowitz fc Johnson 1995 ). One of the most remarkable properties of such layout is that there exist 5 linear combinations of 
the resonators' readouts which are directly proportional to the 5 quadrupole GW amplitudes g^^"^\t) of equation (^). M&J 
call these com binations mode c hannels, and they therefore play a fundamental role in GW signal deconvolution in a real, 
noisy system ( Merkowitz 1998 ; Merkowitz, Lobo fc Serrano 199£ ). In addition, a reduced scale prototype antenna — called 



TIGA, for Truncated Icos ahedron Gravitat ional Antenna — was constructed at Lousiana State University, an d its working 
exper imentall y put to test (Merkowitz 1995). The remarkable success of this experiment in almost every detail ( Merkowitz & 



Johns an 199t ; 



Merkowitz fc Johnson 1997 



Merkowitz fc Johnson 1998 ) stands as a vivid proof of the practical feasibility of 



a spherical GW detector ( Astone et al. 19971 ) . 

Despite its success, the theoretical model proposed by MfcJ to describe the system dynamic s is based upon a simplifying 
assumption that the resonat ors only couple to to the quadrupole vibration modes of the sphere ( Johnson fc Merkowitz 199^ ; 



Merkowitz fc Johnson 1995). While this is see n a posteriori of experimental measurements to be a very good approximation 



( Merkowitz 1995 ; Merkowitz fc Johnson 1997 ), a deeper physical reason which explains why this happens is missing so far. 
The original motivation for the research we present in this article was to develop a new and more gen eral approa ch for the 
analysis of the resonator problem, very much in the spirit of the methodology and results of reference ( Lobo 1995 ); this, we 
thought, would not only provide the necessary tools for a rigorous analysis of the system dynamics, but also contribute to 
improve our understanding of the physics of the spherical GW detector. 

Pursuing this programme, we succeeded in setting up and solving the equations of motion for the coupled system of sphere 
plus resonators. The most important characteristic of our solution is that it is expressible as a perturbative series expansion 
m ascending powers of the small parameter 77^''^, where 77 is the ratio between the average resonator's mass and the sphere's 
mass. The dominant (lowest ) order terms in this expansion appear to exactly reproduce Merkowitz and Johnson's equations 
( Merkowitz fc Johnson 1995 ) , whence a quantitative assessment of their degree of accuracy, as well as of the range of validity 



of their underlying hypotheses obtains; if further precision is required then a well defined procedure of going to next (higher) 
order terms is unambiguously prescribed by the system equations. 

Beyond this, though, the simple and elegant algebra which emerges out of the general scheme has enabled us to explore 



different resonator layo uts, alternative to the unique TIGA of MfcJ. In particular we have found one ( Lobo fc Serrano 1996 



Lobo fc Serrano 1997 ) requiring 5 rather than 6 resonators per quadrupole mode sensed and possessing the remarkable 
property that mode channels can be constructed from the system readouts, i.e., five linear combinations of the latter which 
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Figure 1. Schematic diagramme of the coupling model between a solid sphere and a resonator. The notation is that in the text, but 
subindices have been dropped for clarity. The dashed-dotted arc line on the left indicates the position of the undeformed sphere's surface, 
and the solid arc its actual position. 



are directly proportional to the five quadrupole GW amplitudes. We have called this distribution PHC — see below for full 
details. 

The intrinsically perturbative nature of our approach makes it also particularly well adapted to assess the consequences of 
small defects in the system structure, such as for example symmetry breaking due to suspension attachments, small resonator 
mistunings and mislocations, etc. We have applied them with outstanding success to account for the reported frequency 
measurements of the LSU TIG A prototype ( Merkowitz 199fj ), which was diametrically drilled for suspension purposes; in 
particular, discrepancies between measured and calculated values (generally afi'ecting only the fourth decimal place) are 
precisely of the theoretically predicted order of magnitude. 

We have also applied our methods to analyse the stability of the spherical detector to several mistuned parameters, with 
the result that it is not very sensitive to small construction errors. This conforms again to experimental reports (Merkowitz 
& Johnson 1998), but has the advantage that the argument depends on analytic mathematic al work rather than on compu ter 
simulated procedures, the only ones available to this date to our knowledge — see e.g. (Merkowitz & Johnson 1998) or 
( ^tevenson 1997| ). 

The paper is structured as follows: in section ^ we present the main physical hypotheses of the model, and the general 
equations of motion. In section ^ we set up a Green function approach to solve those equations, and in section ^ we apply it to 
assess the system response to both monopole and quadrupole GW signals. In section |5| we describe in detail the PHC layout, 
including its frequency spectrum and mode channels. Section Jq contains a few brief considerations on the system response 
to a hammer stroke calibration signal, and finally in section m we study how the different parameter mistunings affect the 
detector's behaviour. The paper closes with a summary of conclusions, and three appendices where the heavier mathematical 
details are made precise for the interested reader. 



2 GENERAL EQUATIONS 



With minor improvements, we shall use the notation of references (Lobo 1991:) and (Lobo & Serrano 1996), some of which 



is now briefly recalled. We consider a solid sphere of mass A4, radius R, (uniform) density g, and elastic Lame coefficients A 
and /i, endowed with a set of J resonators of masses Ma and resonance frequencies Qa {a — l,...,J), respectively. We shall 
model the latter as point masses attached to one end of a linear spring, whose other end is rigidly linked to the sphere at 
locations Xa — see Figure The system degrees of freedom are given by the field of elastic displacements u(x, t) of the sphere 
plus the discrete set of resonator spring deformations Za{t); equations of motion need to be written down for them, of course, 
and this is our next concern in this section. 

We shall assume that the resonators only move radially, and also that Classical Elasticity theory (Landau & Lifschitz 
197C) is sufficiently accurate for our purposes^. In these circumstances we have (Lobo & Serrano 1996) 



t We clearly do not expect relativistic motions in extremely small displacements at typical frequencies in the range of 1 kHz. 
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= MV'u+(A + M)V(V-u)+f(x,t) (2) 

■iait) = -'nl[Za{t)~u4t)]+C~\t) , a = l,...,J (3) 

where iia =Xa/R is the outward pointing normal at the the a-th resonator's attachment point, and 

Ua{t) = na-u{xa,t) , a = 1, . . . , J (4) 

is the radial deformation of the sphere's surface at x^. A dot C) is an abbreviation for time derivative. The term in square 
brackets in (^) is thus the spring deformation — q{t) in Figure g. 

f(x, t) in the rhs of (^) contains the density of all non-mternal forces acting on the sphere, which we expediently split 
into a component due the resonators' back action and an external action proper, which can be a GW signal, a calibration 
signal, etc. Thus 

f fx, — frcsonators (x, -j- fcxtcrnal (x, (5) 

Finally, in the rhs of (||) is the force per unit mass (acceleration) acting on the a-th resonator due to external 

agents. 

Since we are making the hypothesis that the resonators are point masses the following holds: 



J 

frcsonators(x, t) = ^ MaO,l [Za{t) - Ua{t)] (5*^^ (x - Xa) Ua (6) 

a = l 

where 6^^^ is the three dimensional Dirac density function. 

The external forces we shall be considering in this paper will be gravitational wave signals, and also a simple calibration 
signal, a perpendicular hammer stroke. GW driving terms, we recall from (^), can be written 

fGw(x,t)=f(0°)(x)5(°«)(t) + ^ f(^'">(x)p(2'">(t) (7) 

m—~2 



for a general metric wave — see (Lobo 1995) for explicit formulas and technical details. While the spatial coefficients f (x) 
are pure form factors associated to the tidal character of a GW excitation, it is the time dependent factors g^''"^\t) which 
carry the specific information on the incoming GW. The purpose of a GW detector is to determine the latter coefRcients on 
the basis of suitable measurements. 

If a GW sweeps the observatory then the resonators themselves will also be affected, of course. They will be driven, 
relative to the sphere's centre, by a tidal acceleration which, since they only move radially, is given by 

= R0i0j{t) Xa,ina,j , a = 1, . . . , J (8) 

where Roiojit) are the "electric" components of the GW Riemann tensor at the centre of the sphere. These can be easily 
manipulated to giv^ 

yUr^a)g^'^\t) , a = l,...,J (9) 

l — O and 2 
m — — Z , . . . ,Z 

where 7? is the sphere's radius. 

We shall also be later considering in this paper the response of the system to a particular calibration signal, consisting in 
a hammer stroke with intensity fo, delivered perpendicularly to the sphere's surface at point xq: 

fstrokc(x,t) =fo5'=^>(x-Xo)5(i) (10) 

which we have modeled as an impulsive force in both space and time variables. Unlike GW tides, a hammer stroke will be 
applied on the sphere's surface, so it has no direct effect on the resonators. In other words, 

Cr°'^(i)=0, a=l,...,J (11) 
Our fundamental equations thus finally read: 



■t Yjm(n) are spherical harmonics (Edmonds 1960) — see also the multipole expansion of RoiOjit) in reference (Lobo 1995). 
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d u 

= A'V'u+(A + fi)V(V-u) + 
J 

- Xi,) lit, + fcxternal 

(x,t) (12) 

6=1 

■ia{t) = -nl[Zait)-Ua{t)]+C~\t) , a=l,...,J (13) 

where fcxtcrnai(x, i) will be given by either (Q) or (p^, as the case may be. Likewise, ^™*°'^"^'(t) will be given by (^) or 
(^l|), respectively. The remainder of this paper will be concerned with finding solutions to the system of coupled differential 
equations ( |l^ ) and (|l^), and with their meaning and consequences. 



3 GREEN FUNCTION FORMALISM 



In order to solve equations (p^-(|l^) we shall resort to Green function formalism. The essentials of this procedure in the context 
of the present problem can bo found in detail in reference (Lobo 1995); more specific technicalities are given in appendix ^ 



By means of such formalism equations (p^-(p^) become the following integro-differential system 



J 



Ua{t) = uT"''''^\t)+J2r!t Kab(t-t') [zt(t')-ub{t')] dt' (14) 

6=1 "'O 



Ut) = C (t)-f^a , a = l,...,J (15) 

where ul^^°'""^\t) = na ■ u'"''°™'''(xa, t), and u'""°™'''(x, t) is the bare (i.e., without attached resonators) sphere's response to 
the external forces foxtcrnai(x, t) in the rhs of (p^). Kab{t) is a kernel matrix defined by the following weighted sum of diadic 
products of wavefunctions^ 

Kabit) = ^ i^N^ [ni,-u^(x6)] [na-ujv(xa)] sinwjvt (16) 

N 

Finally, we have defined the mass ratios of the resonators to the entire sphere 

^i>=^, 6=1,..., J (17) 

which will be small parameters in a real device. 

Before proceeding further, let us briefly pause for a qualitative inspection of equations (|l|) and (Q. Equation ^ 
shows that the sphere's surface deformations Ua{t) are made up of two contributions: one due to the action of external agents 
(GWs or other), contained in ■it™*'^"^"'''(t), and another one due to coupling to the resonators. The latter is commanded by the 
small parameters rjb, and correlates to all of the sphere's spheroidal eigenmodes through the kernel matrix Kabit). This has 
consequences for GW detectors, for even though GWs may only couple to quadrupole and monopolep[] spheroidal modes of 
the free sphere (Lobo 1995; Bianchi et al. 199^ , attachment of resonators causes, as we see, coupling between these and the 



other modes of the antenna, and conversely, these modes back-act on the former. As we shall shortly prove, such effects can 
be minimised by suitable tuning of the resonators' frequencies. 



3.1 Laplace transform domain equations 

We now take up the problem of solving these equations. Equation (^) is an integral equation belonging in the general 
category of Volterra equations (Tricomi 1957), but the usual iterative solution to it by repeated substitution of Ubit) into 



the kernel integral is not viable here due to the dynamical contribution of Zb(t), which is in turn governed by the 
equation (|l^). 

A better suited method to solve this integro-differential system is to Laplace-transform it. We denote the Laplace transform 
of a generic function of time f(t) with a caret [") on its symbol, e.g.. 



3 We shall often use the capitalised index N to imply the multiple index {nlm} which chara cterises the sphere's w avcfunctions. 

^ Monopole modes only exist in scalar-tensor theories of gravity, such as e.g. Brans— Dicke (Brans & Dicke 1961); General Relativity 

does not belong in this category. 
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f{s)= / f{t)e-''dt (18) 
Jo 

and make the assumption that the system is at rest before an instant of time, t = 0, say, or 

U(x,0)=u(x,0)=2a(0)=ia(0)=0 (19) 

Equations ( |l^ ) and ( [l5| ) then adopt the equivalent form 
J 

Ua{s) = Mr*''""'"'(a) - y^^VbKabis) [zbjs) - Ubjs)] (20) 
b=l 

s'^ Za{s) = ^^^""'''(s) - f^a [5a(s) - tia(s)] , a = l,...,J (21) 

for which use has been made of tlie convolution theorem for Laplace transforms|[J]. A further simplification is accomplished if 
we consider that we shall in practice be only concerned with the measurable quantities 

qa{t) = Za{t) - Ua{t) , a=l,...,J (22) 

representing the resonators' actual elastic deformations — cf. Figure 0. It is readily seen that these verify the following: 



E 



Sab + Vb ^2 _^ ^2 Kab{s) 



2 (^external / \ 

" / \ '-' "external / \ , Sa i t /r»o\ 

g''^^) ^" 524.^2 "a {S)+ ^2j^^2 ' a=l,...,J (23) 



Equations (^) constitute a significant simplification of the original problem, as they are a set of just J algebraic rather 
than integral or differential equations. We must solve them for the unknowns qa{s), then perform inverse Laplace transforms 
to revert to qa{t). We do this next. 



4 SYSTEM RESPONSE TO A GRAVITATIONAL WAVE 

Our concern now is the actual system response when it is acted upon by an incoming GW. We shall calculate it by making a 
number of simplifying assumptions, more precisely: 

i) The detector is perfectly spherical. 

ii) The resonators have identical masses and resonance frequencies. 

iii) The resonators' frequency is accurately matched to one of the sphere's oscillation eigenfrequencies. 

As we shall see below (section ^), a real system can be appropriately treated as one which deviates by definite amounts 
from this idealised construct. Therefore detailed knowledge of the ideal system behaviour is essential for all purposes. This is 
the justification for the above simplifications. 

The wavefunctions u„im,(x) of an elastic sphere can be found in reference (Lobo 1995) in full detail, and we shall keep 
the notation of that paper for them. The Laplace transform of the kernel matrix (|l6[) can "thus be expressed as — see equation 
(A21) in appendix^: 

^-(^)=E ^^'("-"^)-E i^xir (24) 

nl nl 

where the last term simply defines the quantities xib''- Note that the sums here stretch across the entire spectrum of the 
solid sphere. 

Our next assumption that all the resonators are identical simply means that 

rji = . . . = r]j = rj , D,i = . . . = D,j = D, (25) 

The third hypothesis makes reference to the fundamental idea behind using resonators, which is to have them tuned to 
one of the frequencies of the sphere's spectrum. We express this by 

^ = w„oio (26) 



I This theorem states, it is recalled, that the Laplace transform of the convolution product of two functions is the arithmetic product 
of their respective Laplace transforms. 
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where LiJnoiQ is a specific and fixed frequency of tfie spheroidal spectrum. 

In a GW detector it will only make sense to choose lo = or lo — 2, as only monopole and quadrupole sphere modes 



couple to the incoming signal; in practice, no will refer to the first or perhaps second harmonic (Coccia et al. 1995a). We shall 
however keep the generic expression ( p^ for the time being in order to encompass all the possibilities with a unified notation. 
Based on the above hypotheses, equation (ESh can be rewritten in the form 



J 

E 

6=1 



Sab + ri 



nl 



=- ^24,5^2 ^° (g)+ s2 + Q2 ' i^^^nolo) (27) 



where £,a (s) is the Laplace transform of (||), i.e., 

i2^is)^R ^ F,„K)p("")(s) , a = l,...,J (28) 

l — O and 2 
in — — l,...,l 

As mentioned at the end of the previous section, we must now invert the matrix in the Ihs of (^^, which will give us an 
expression for qa{s), then find the inverse Laplace transform of these functions to revert back to the time domain. A simple 
glance at the equation suffices however to grasp the unsurmountable difficulties of accomplishing this analytically. 

Thankfully, though, a perturbative approach is applicable when the masses of the resonators are small compared to the 
mass of the whole sphere, i.e., when the inequality 

7? < 1 (29) 

holds. We shall henceforth assume that this is the case, as also is with cylindrical bar resonant transducers. It is shown in 
appendix ^ that the perturbative series happens in ascending powers of r]^^^, rather than rj itself, and that the lowest order 
contribution has the form 

ga(s) = '7"'^' 51 + ^(°) ' <^ = h---,J (30) 

l,m 

where 0(0) stands for terms of order rj" or smaller. Here, Ai' (s; Q,) is a transfer function matrix which relates linearly the 
system response qa{s) to the GW amplitudes g''"''(s), in the usual sense that qa{t) is given by the convolution product of the 
signal (f) with the time domain expression, Ai''"'(t;n), of Ai'™'(s;n). The detector is thus seen to act as a linear filter 
on the GW signal, whose frequency response is characterised by the properties of K^a"^\s; Q,). More specifically, the filter has 
a number of characteristic frequencies which correspond to the imaginary parts of the poles of Aa'"'(s; f2). As also shown in 
appendix pi these frequencies are the symmetric pairs 



2 



(^l±y^^ lA„o,„(7?)|Ca7)^/'j +0(?7) , a = l,...,J (31) 
where C,a is the a-th eigenvalue of the Legendre matrix 

Pig{na-nb) , a, 6=1, ...,J (32) 

associated to the multipole {Iq) selected for tuning — see (p6|). These frequency pairs correspond to heats, typical of resonantly 
coupled oscillating systems — we shall find them again in section ^ in a particularly illuminating example. 

Another very important fact is also neatly displayed by equation (^^: the resonators' motions are mechanically amplified 
by a factor ri~^^'^ relative to the driving amplitudes 3''™'(s). This is the counterpart, in our multimode system, of a similar 
behaviour known to happen in monomode cylindrical antennas ( Astone et al. 19'9^ ). 



The specific form of the transfer function matrix Aa'"-'(s;f2) depends on both the selected mode to tune the resonator 
frequency Q., and on the resonator distribution geometry. We now come to a discussion of these. 



4.1 Monopole gravitational radiation sensing 



General Relativity, as is well known, forbids monopole GW radiation. More general metric theories, e.g. Brans-Dicke (Brans & 



Dicke 1961), do however predict this kind of radiation. It appears that a spherical antenna is potentially sensitive to monopole 
waves, so it can serve the purpose of thresholding, or eventually detecting them. This clearly requires that the resonator set 
be tuned to a monopole harmonic of the sphere, i.e., 

Q = uj^o , {lo = 0) (33) 

where n tags the chosen harmonic — most likely the first {n = 1) in a thinkable device. 
Since Po{z) = 1 (for all 2) the eigenvalues of Po(na-n6) are, clearly. 
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J , C2 = • • . = C5 = 

for any resonator distribution. The tuned mode frequency thus sphts into a single strongly coupled pair: 



(34) 



2 



1 ± 



\A„,{R)\ri'^'\+Oiv) 



(35) 



The A-matrix of equation 
K '{s;ujno) = 
whence the system response is 



is seen to be in this case 



f 2 , 2 



r-i 



2 , 2 - 

S + LU_ 



-1/2 



(-1)^ 



5 in Sr, 



1,. 



(36) 



(37) 



regardless of resonator positions. The overlap coefficient a„o is calculated by means of formulas given in (Lobo 1995), and has 
dimensions of length. By way of example, aio/-R = 0.214, and 020/ R~ —0.038 for the first two harmonics. 

A few interesting facts are displayed by equation (|^. First, as we have already stressed, it is seen that if the resonators 
are tuned to a monopole detector frequency then only monopole wave amplitudes couple strongly to the system, even if 
quadrupole radiation amplitudes are significantly high at the observation frequencies uj±. Also, the amplitudes ga(s) are 
equal for all a, as corresponds to the spherical symmetry of monopole sphere's oscillations, and are proportional to J~^^^, a 
factor we should indeed expect as an indication that GW energy is evenly distributed amongst all the resonators. A single 
transducer suffices to experimentally determine the only monopole GW amplitude g'°°-'(s), of course, but (^^ provides the 
system response if more than one sensor is mounted on the antenna for whatever reasons. 



4.2 Quadrupole gravitational radiation sensing 

We now consider the more interesting case of quadrupole motion sensing. We thus take 

Q = Ldn2 , (^0 = 2) 



(38) 



where n labels the chosen harmonic — most likely the first (n = 1) or the second (n = 2) in a practical system. The evaluation 
of the A-matrix is now considerably more involved (Serrano 1999), yet a remarkably elegant form is found for it: 



Ai-)(s;c.„.) = (-l)-Jf a„.5] K] 



2, 2\-l ('2, 2\-l 
S - [S +UJc-) 



(c) (c)* 



y2m(n6) Si: 



(39) 



where vi'^^ is the c-th normalised eigenvector of P2(na-n6), associated to the non-null eigenvalue Let us stress that equation 
( pgj ) explicitly shows that at most 5 pairs of modes, of frequencies uJc±, couple strongly to quadrupole GW amplitudes, no 
matter how many resonators in excess of 5 are mounted on the sphere. The tidal overlap coefficients a2n can also be calculated, 
and give for the first two harmonics (Lobo & Serrano 1996) 



ai2 



0.328 



a22 
R 



0.106 



(40) 



The system response is thus 



''=1 LCcT^O 



(c) (c)* 



J2 Y^rnin,) g(^^\s) + O{0) , a = l...,J 



(41) 



Equation (|4l[ ) is completely general, i.e., it is valid for any resonator configuration over the sphere's surface, and for any 
number of resonators. It describes precisely how all 5 GW amplitudes g'^'"'(s) interact with all 5 strongly coupled system 
modes; like before, only quadrupole wave amplitudes are seen in the detector (to leading order) when Q = lj„2 , even if the 
incoming wave carries significant monopole energy at the frequencies lOc±. 

The degree of generality and algebraic simplicity of (^) is new in the literature. As we shall now see, it makes possible 
a systematic search for different resonator distributions and their properties. 
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5 THE PHC CONFIGURATION 



Merkowitz and Johnson's TIG A (Johnson & Merkowitz 1993) is highly symmetric, and is the minimal set with maximum 
degeneracy, i.e., all the non-null eigenvalues C^a are equal. To accomplish this, however, 6 rather than 5 resonators are required 
on the sphere's surface. Since there are just 5 quadrupole GW amplitudes one may wonder whether there are alternative 
layouts with only 5 resonators. Equation ( pl| ) is completely general, so it can be searched for an answer to this question. In 
reference ( |Lobo fc Serrano 1996 ) we made a specific proposal, which we now describe in more detail. 

In pursuing a search for 5 resonator sets we found that distributions having a sphere diameter as an axis of pentagonal 
symmetry^^ exhibit a rather appealing structure. More specifically, let the resonators be located at the spherical positions 



(all a) 



VPa = (a - 1) — 
5 



a = 1, 



(42) 



The eigenvalues and eigenvectors of P2(na-ni,) are easily calculated: 



Co' = I (3 cos^ Q - 1) 



(m) 



■ C,Tr^ Y-im{ja.a) , m = — 2, ...,2, a = l 
so the A-matrix is also considerably simple in structure in this case: 



yi _ yi _ 15 
',2 - t,-2 - — 

5 



47r 1 



I 2 , 2\-l ( 1 , 2\-l 



Ai'™)(s;a;„2) = 
where we have used the notation 

= n' \Ar.2{R)\ C™ n^'^ j + 0(77) , m = -2, 



y2m(n,)5,2 (pnc) 



(43) 
(44) 

(45) 
(46) 



As we see from these formulas, the five expected pairs of frequencies actually reduce to three, so pentagonal distributions 
keep a certain degree of degeneracy, too. The most important distinguishing characteristic of the general pentagonal layout is 
best displayed by the explicit system response: 



Qa(s) = -77" 



-1/2 /IF 



an2 



+ 



+ 



1 

2Co 
1 

20 
1 

%2 



y2o(n<.)5('°>(s) 



{s'+Ji^) '-(s' '] [y2i(n„)p("'(s) + r2-i(n.)5(i-^)(s)] 
(s'+c.i+)"'- (s'+^2'-)"'l [F22(n„)5*''^(s) + y2-2(n„)p('-')(s)] 



(47) 



This equation indicates that different wave amplitudes selectively couple to different detector frequencies. This should be 
considered a very remarkable fact, for it thence follows that simple inspection of the system readout specfrum^ immediately 
reveals whether a given wave amplitude (f"^{s) is present in the incoming signal or not. 

Pentagonal configurations also admit mode channels, which are easily constructed from (^) thanks to the orthonormality 
property of the eigenvectors (^): 



" (m) / \ 



E(m)* - 

a = l 



qa{s) = r) a„2 



/ 2 , 2 



( 2 , 2 



g^^^\s) + 0{Q) 



m = -2, 



(48) 



These are almost identical to the TIG A mode channels ( Merkowitz fc Johnson 1995 ) , the only difference being that each 
mode channel comes now at a single specific frequency pair ujm±- 



Mode cha nnels are fundamental in signal deconvolution algorithms in noisy systems ( Merkowitz 1998 ; Merkowitz, Lobo fc 
Serrano 1999). Pentagonal resonator configurations should thus be considered non-trivial candidates for a real GW detector. 

Based on these facts one may next ask which is a suitable transducer distribution with an axis of pentagonal symmetry. In 
Figure ^ we give a plot of the eigenvalues (^sj) as a function of a, the angular distance of the resonator set from the symmetry 
axis. Several criteria may be adopted to select a specific choice in view of this graph. An interesting one was proposed by 
us in reference (Lobo & Serrano 1996) with the following argument. If for ease of mounting, stability, etc., it is desirable to 



By this we mean resonators are placed along a parallel of the sphere every 72° . 
tt In a noiseless system, of course 
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Figure 2. The three distinct eigenvalues Cm (r>i = 0.1.2) as functions of the distance of the resonator parallel's co-latitude a relative to 
the axis of symmetry of the distribution, cf. equation (flSh. 





Figure 3. To the left, the pentagonal hexacontahedron shape. Certain faces are marked to indicate resonator positions in a specific 
proposal — see text — as follows: a square for resonators tuned to the first quadrupole frequency, a triangle for the second, and a star 
for the monopole. On the right we see the (pentagonal) face of the polyhedron. A few details about it: the confluence point of the dotted 
lines at the centre is the tangency point of the inscribed sphere to the PHC; the labeled angles have values = 61.863°, /3 = 87.205°; the 
angles at the T-vertices are all equal, and their value is 118.1366°, while the angle at P is 67.4536°; the ratio of a long edge (e.g. PTi) 
to a short one (e.g. T\T2) is 1.74985, and the radius of the inscribed sphere is twice the long edge of the pentagon, R = 2 PT\. 
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Figure 4. Compared line spectrum of a coupled TIGA and a PHC resonator layout in an ideally spherical system. The weakly coupled 
central frequency in the TIGA is drawn dashed. The frequency pair is 5-fold degenerate for this layout, while the two outer pairs of the 
PHC are doubly degenerate each, and the inner pair is non-degenerate. Units in abscissas are r)^/^Q, and the central value, labeled 0.0, 
corresponds to Q. 

have the detector milled into a close-to-spherical polyhedric shapej^ then polyhedra with axes of pentagonal symmetry must 



be searched. The number of quasi regular convex polyhedra is of course finite — there actually are only 18 of them (Holden 



19771 

we 



Har'El 1993) — , and we found a particularly appealing one in the so called pentagonal hexacontahedron {PHC), which 
see in Figure ^ , left. This is a 60 face polyhedron, whose faces are the identical irregular pentagons of Figure ^, right. The 
PHC admits an inscribed sphere which is tangent to each face at the central point marked in the Figure. It is clearly to this 
point that a resonator should be attached so as to simulate an as perfect as possible spherical distribution. 

The PHC is considerably spherical: the ratio of its volume to that of the inscribed sphere is 1.057, which quite favourably 
compares to the value of 1.153 for the ratio of the circumscribed sphere to the TI volume. If we now request that the frequency 
pairs ujm± be as evenly spaced as possible, compatible with the PHC face orientations, then we must choose a = 67.617° , whence 

t^o± =t^i2 (l± 0.5756 r;^/^) , uji± ^ lui2 {l ± 0.8787 rj^^^) , a;2± = t^i2 (l ± 1.0668 j?^/^) (49) 

for instance for n = a;i2, the first quadrupole harmonic. In Figure ^ we display this frequency spectrum together with the 
multiply degenerate TIGA for comparison. 

The criterion leading to the PHC proposal is of course not unique, and alternatives can be considered. For example, 
if the 5 faces of a regular icosahedron are selected for sensor mounting (a = 63.45°) then a four- fold degenerate pair plus a 
single non-degenerate pair is obtained; if the resonator parallel is 50° or 22.6° away from the "north pole" then the three 
frequencies ujo+, and UJ2+ are equally spaced; etc. The number of choices is virtually infinite if the sphere is not milled 



into a polyhedric shape (Stevenson 1997; Magalhaes et al. 1997). 

Let us finally recall that the complete PHC proposal (Lobo & Serrano 1996| ) was made with the idea of building an as 



complete as possible spherical GW antenna, which amounts to making it sensitive at the first two quadrupole frequencies and 
at the first monopole one. This would take advantage of the good sphere cross section at the second quadrupole harmonic 



(Coccia et al. 1995a), and would enable measuring (or thresholding) eventual monopole GW radiation. Now, the system 
pattern matrix Aji"" (s; f2) has identical structure for all the harmonics of a given I series — see ( |36| ) and ( |39| ) — , and so too 
identical criteria for resonator layout design apply to either set of transducers, respectively tuned to uj-12 and 0^22 • The PHC 
proposal is best described graphically in Figure ^, left: a second set of resonators, tuned to the second quadrupole harmonic cj22 
can be placed in an equivalent position in the "southern hemisphere" , and an eleventh resonator tuned to the first monopole 



M This is the philosophy suggested and experimentally implemented by Merkowitz and Johnson at L3U. 
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frequency u>io is added at an arbitrary position. It is not difficult to see, by the general methods outlined earlier on in this 
paper, that cross interaction between these three sets of resonators is only second order in r]^^^, therefore weak. 

A spherical GW detector with such a set of altogether 11 transducers would be a very complete multi-mode multi-frequency 
device with an unprecedented capacity as an individual antenna. Amongst other it would practically enable monitoring of 
coalescing binary chirp signals by means of a rather robust double passage method (Coccia & Fafone 199C), a prospect which 



was considered so far possible only with broadband long baseline laser interferometers (Dhurandhar, Krolak & Lobo 198£ 



Krolak, Lobo & Meers 1991), and is almost unthinkable with currently operating cylindrical bars. 



6 A CALIBRATION SIGNAL: HAMMER STROKE 

This section is a brief digression from the main streamline of the paper. We propose to assess now the system response to a 
particular, but useful, calibration signal: a perpendicular hammer stroke. 

We first go back to equation (ESI) and replace ■it™*'^"'"'''(s) in its rhs with that corresponding to a hammer stroke, which 
is easily calculated — cf. appendix ^ 

fo 



"Stroke/ \ 



nl 



\A„i{R)\^ Pi{na-no) , a = 1, . . . , J 



(50) 



where no are the spherical coordinates of the hit point on the sphere, and fo = no-fo/A1. Clearly, the hammer stroke excites 
all of the sphere's vibration eigenmodes, as it has a completly flat spectrum. 

The coupled system resonances are again those calculated in appendix The same procedures described in section ^ for 
a GW excitation can now be pursued to obtain 



Us) = v-'^'i-i) 



21 + 1 
4% 



fo \An,{R)\ X 



E EM(^' 



6=1 kCcT^O 



(c) (c)* 



Pi{nfno) + OiO) 



a = 1 



(51) 



when the system is tuned to the nl-th spheroidal harmonic, i.e., n — cuni- It is immediately seen from here that the system 
response to this signal when the resonators are tuned to a monopole frequency is given by 



q4s) = ^-'^^ (-1) 



7^ \A,MR)\ i 



2 \-l 



(52) 



an expression which holds for all a, and is independent of either the resonator layout or the hit point, which in particular 
prevents any determination of the latter, as obviously expected. The frequencies uj± are those of (^), and we find here again 
a global factor J~^^^, as also expected. 

We consider next the situation when quadrupole tuning is implemented, Q = uj„2. We shall however do so only for the 
PHC and TIGA configurations, as more general considerations are not quite as interesting at this point. 

6.1 PHC and TIGA response to a hammer stroke 

Expanding equation ^ by substitution of the eigenvalues (rn and eigenvectors iii™^ of the PHC, one readily finds that the 
system response is given by 



qa{s) = V-'^\foJ^\A„2{R)\ J2 \ 



( 2 , 2 



-1 ( 2 , 2 \-l 



Y2 (na) Y2 

m 

(no) , 



a — 1, 



(53) 



and the mode channels by 

y('"'(s) =77-^/^0 I^n2(i?)! 



>^2m(no) , m = -2, 



(54) 



These equations indicate that the system response qa{t) is a superposition of three different beat^^ while the mode 
channels are single beats each, but with differing modulation frequencies. This is represented graphically in Figure |5|, where 



§§ A beat is a modulated oscillation of the form sin i (a;+ — i^}- )t cos Ht, where a;+ and ijJ— are nearby frequencies, and c<;+ + u)- = 2 f2. 



The Laplace transform of such function of time is precisely (fi/2) 
difference aj+ —ijJ-, which in our case is proportional to r;^/^. 



up to higher order terms in the 



© 1999 RAS, MNRAS 000, [^-^ 




Figure 5. Simulated response of a PHC to a liammer stroke: the time series and their respective spectra, both for direct resonator 
readouts and mode channels. Note that while the former are not simple beats, the latter are. 
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we see the result of a numerical simulation of the PHC response to a hammer stroke, delivered to the solid at a given location. 
The readouts qa{t) are somewhat complex time series, whose frequency spectrum shows three pairs of peaks — in fact, the 
lines in the ideal spectrum of Figure ^ The mode channels on the other hand are pure beats, whose spectra consist of the 
individually separate pairs of the just mentioned peaks. 

The response of the TIGA layout to a hammer stroke has been described in detail by Merkowitz and Johnson — see e.g. 



reference (Merkowitz & Johnson 1997). Our formalism does of course recover the results obtained by those authors; in the 



notation of this paper, we have 



(la(S) 



y'- '{s) 



,-1/2 



/|^/olA.(i?)|i 



( 2 , 2 \ -1 
S +UJ+) 



12, 2 \ -1 



P2(na-no) 



"i1>m(no) 



TIGA 
,2 



(55) 
(56) 



for the system response and the mode channels, respectively, where 



2 2 
^± = ^n2 



1±\/^ \An2{R)W'^ 



oiv) 



(57) 



are the five-fold degenerate frequency pairs corresponding to the TIGA distribution. Comparison of the mode channels shows 
that they are identical for PHC and TIGA, except that the former come at different frequencies depending on the index m. 
One might perhaps say that the PHC gives rise to a sort of "Zeeman splitting" of the TIGA degenerate frequencies, which 
can be attributed to an axial symmetry breaking of that resonator distribution: the PHC mode channels partly split up the 
otherwise degenerate multiplet into its components. 



7 SYMMETRY DEFECTS 

So far we have made the assumption that the sphere is perfectly symmetric, that the resonators are identical, that their 
locations on the sphere's surface are ideally accurate, etc. This is of course unrealistic. So we propose to address now how 
such departures from ideality affect the system behaviour. As we shall see, the system is rather robust, in a sense to be made 
precise shortly, against a number of small defects. 

In order to quantitatively assess ideality failures we shall adopt a philosophy which is naturally suggested by the results 
already obtained in an ideal system. It is as follows. 

As we have seen in previous sections, the solution to the general equations ( psj ) must be given as a perturbative series 
expansion in ascending powers of the small quantity rj^^^. This is clearly a fact not related to the system's symmetries, so it 
will survive symmetry breakings. It is therefore appropriate to parametrize deviations from ideality in terms of suitable powers 
of rj^^^, in order to address them consistently with the order of accuracy of the series solution to the equations of motion. An 
example will better illustrate the situation. 

In a perfectly ideal spherical detector the system frequencies are given by equations (^l|). Now, if a small departure from 
e.g. spherical symmetry is present in the system then we expect that a correspondingly small correction to those equations 
will be required. Which specific correction to the formula will actually happen can be qualitatively assessed by a consistency 
argument: if symmetry defects are of order rj^^^ then equations (^l|) will be significantly altered in their r/^^^ terms; if on the 
other hand such defects are of order 77 or smaller then any modifications to equations (^) will be swallowed into the 0(0) 
terms, and the more important r}^^^ terms will remain unaff'ected by the symmetry failure. We will say in the latter case that 
the system is robust to that ideality breaking. 

More generally, this argument can be extended to see that the only system defects standing a chance to have any influences 
on lowest order ideal system behaviour are defects of order r}^^^ relative to an ideal configuration. Defects of such order are 
however not necessarily guaranteed to be significant, and a specific analysis is required for each specific parameter in order to 
see whether or not the system response is robust against the considered parameter deviations. 

We therefore proceed as follows. Let P be one of the system parameters, e.g. a sphere frequency, or a resonator mass or 
location, etc. Let Pidcai be the numerical value this parameter has in an ideal detector, and let Prcai be its value in the real 
case. These two will be assumed to differ by terms of order r/^^^, or 

Prcal = Pidcai (I+P77'/') (58) 

For a given system, p is readily determined adopting (^^ as the definition of Proai, once a suitable hypothesis has been 
made as to which is the value of Pidcai- In order for the following procedure to make sensible sense it is clearly required that 
p be of order 1 or, at least, appreciably larger than rj^^'^ ■ Should p thus calculated from (js^ ) happen to be too small, i.e., of 
order 77^^^ itself or smaller, then the system will be considered robust as regards the affected parameter. 

We now apply this criterion to various departures from ideality. 
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7.1 The suspended sphere 



An earth based observatory obviously requires a suspension mechanism for the large sphere. If a nodal point suspension is 
e.g. selected then a diametral bore has to be drilled across the sphere (Merkowitz 1995). The most immediate consequence of 
this is that spherical symmetry is broken, what in turn results in degeneracy lifting of the free spectral frequencies u!„i, which 
now split up into multiplets uj„im {m — —l,...,l). The resonators' frequency Q cannot therefore be matched to the frequency 
, K,,f at mr^ot +^ ojje of the Ungigrn's. lu this subsection we keep the hypothesis that all the resonators are identical — we 
and assume that falls within the span of the multiplet of the cOnoiom's. Then we write 

l/2x 



shall relax it later- 



m = ~/o, 



, ^0 



(59) 



We now search for the coupled frequencies, i.e., the roots of equation (BE). The kernel matrix Kab{s) is however no longer 
given by (^^, due the removed degeneracy of and we must stick to its general expression (A18), or 



^ S2 

nlm 



9/ I 1 

47r 



' nlm 7 ' nlm 

nlm 

Following the steps of appendix ^ we now seek the roots of the equation 

c\1 2 2 



det 



(nlm) 
Xab 



Sab + V 



lo 

E 

m— — Zi 



{s^ + n'^){s^ + ^l^,^j 



(nolom) 
Xab "t" '/ 



nl^nQiQ ,m 



{nlm) 
Xab 



(60) 



(61) 



Since Q. relates to Un^m through equation ( [59| ) we see that the roots of (^l|) fall again into either of the two categories 
(|B9|)-(B10)(see Appendix]^, i.e., roots close to ±iQ. and roots close to ±iijjnim (nlj^nolo). We shall exclusively concentrate 
on the former now. Direct substitution of the series ( [B9| ) into ( |6l| ) yields the following equation for the coefficient xi- 



det 



Sab 



- E 



Xab 



= 



(62) 



when there is full degeneracy. 

Rather, there are 2Zo + l+J of them, with a maximum 



This is a variation of (Bll), to which it reduces when pm =0, i.e 

The solutions to no longer come in symmetric pairs, like (js^ 
number of 2(2Zo+l) non-identically zero roots if J>2Zo + lP|. For example, if we choose to select the resonators' frequency 
close to a quadrupole multiplet (Zo = 2) then ( |62| ) has at most 5+ J non-null roots, with a maximum ten no matter how many 
resonators in excess of 5 we attach to the sphere. Modes associated to null roots of can be seen to be weakly coupled, just 
like in a free sphere, i.e., their amplitudes are smaller than those of the strongly coupled ones by factors of order rf'^^ . 

In order to assess the reliability of this method we have applied it to see what are its predictions for a real system. To this 
end, data taken with the TIGA prototype at LSI7[|||] were used to confront with. The TIGA was drilled and suspended from 
its centre, so its first quadrupole frequency split up into a multiplet of five frequencies. Their reportedly measured values are 



Wl20 



: 3249 Hz , LJ121 = 3238 Hz 



t^l2- 



3236 Hz 



'i'122 



3224 Hz , UJ12 . 



3223 Hz 



All 6 resonators were equal, and had the following characteristic frequency and mass, respectively; 
1 



n = 3241 Hz 



V = 



1762.45 

Substituting these values into (^) it is seen that 
Po = 0.2075 , pi = -0.0777 , p_i = -0.1036 , pa = -0.4393 , p_ 



-0.4650 



(63) 



(64) 



(65) 



Equation (|62|) can now be readily solved once the resonator positions are fed into the matrices xib^™'- Such positions 



correspond to the pentagonal faces of a truncated icosahedron. Merkowitz (Merkowitz 1995) gives a complete account of all 
the measured system frequencies as resonators are progressively attached to the selected faces, beginning with one and ending 
with six. In Figure ^ we present a graphical display of the experimentally reported frequencies along with those calculated 
theoretically by solving equation (^^. In Table 1 we give the numerical values. As can be seen, coincidence between our 
theoretical predictions and the experimental data is remarkable: the worst error is 0.2%, while for the most part it is below 
0.1%. Thus discrepancies between our theoretical predictions and experiment are of order r], as indeed expected — see (p4j). 



In addition, it is also reported in reference (Merkowitz & Johnson 1997) that the 11-th, weakly coupled mode of the TIGA 



This is a mathematical fact, whose proof is relatively cumbersome, and will be omitted here; we just mention that it has its origin 
in the linear dependence of more than 2Zo-|-l spherical harmonics of order Iq. 

Ill These data are contained in reference (Merkowitz 1995), and we want to express our gratitude to Stephen Merkowitz for kindly 
handing them to us. 
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Table 1. Numerical values of measured and theoretically predicted frequencies (in Hz) for the TIGA prototype with varying number of 
resonators. Percent differences are also shown. The calculated values of the tuning and free multiplet frequencies are taken by definition 
equal to the measured ones, and quoted in brackets. In square brackets the frequency of the weakly coupled sixth mode in the full, 
6 resonator TIGA layout. These data are plotted in Figure H. 

Item Measured Calculated % difference Item Measured Calculated % difference 



Tuning 3241 (3241) (0.00) 

Free multiplet 3223 (3223) (0.00) 

3224 (3224) (0.00) 

3236 (3236) (0.00) 

3238 (3238) (0.00) 

3249 (3249) (0.00) 

1 resonator 3167 3164 -0.08 

3223 3223 0.00 

3236 3235 -0.02 

3238 3237 -0.02 

3245 3245 0.00 

3305 3307 0.06 

2 resonators 3160 3156 -0.13 

3177 3175 -0.07 

3233 3233 0.00 

3236 3236 0.00 

3240 3240 0.00 

3302 3303 0.03 

3311 3311 0.00 

3 resonators 3160 3155 —0.15 

3160 3156 -0.13 

3191 3190 -0.02 

3236 3235 -0.02 

3236 3236 0.00 

3297 3299 0.08 

3310 3311 0.02 

3311 3311 0.00 



4 resonators 3159 3155 -0.12 

3160 3156 -0.11 

3168 3165 -0.12 
3199 3198 -0.05 
3236 3236 0.00 
3285 3286 0.03 

3310 3310 0.00 

3311 3311 0.00 
3319 3319 0.00 

5 resonators 3152 3154 0.08 

3160 3156 -0.14 

3163 3162 -0.03 

3169 3167 -0.08 
3209 3208 -0.02 
3268 3271 0.10 
3304 3310 0.17 
3310 3311 0.03 
3313 3316 0.10 
3319 3321 0.06 

6 resonators 3151 3154 0.11 

3156 3155 -0.03 

3162 3162 0.00 

3167 3162 -0.14 

3170 3168 -0.07 
[3239] [3241] [0.06] 
3302 3309 0.23 
3308 3310 0.06 

3312 3316 0.12 
3316 3317 0.02 
3319 3322 0.10 



(enclosed in square brackets in Table 1) has a practically zero amplitude, again in excellent agreement with our general 
theoretical predictions about modes beyond the tenth — see paragraph after equation (js^). 

This is an encouraging result which motivated us to try a better fit by estimating the next order corrections, i.e., xi of 
(B9|). As it turned out, however, matching between theory and experiment did not consistently improve. This is not really that 



surprising, though, as M&J explicitly state (Merkowitz & Johnson 1997) that control of the general experimental conditions in 



which data were obtained had a certain degree of tolerance, and they actually show satisfaction that ~1% coincidence between 
theory and measurement is comfortably accomplished. But 1% is two orders of magnitude larger than rj — cf. equation (|^ — , 
so failure to refine our frequency estimates to order r] is again fully consistent with the accuracy of available real data. 



A word on a technical issue is in order. Merkowitz and Johnson's equations for the TIGA (Johnson & Merkowitz 1995 



Merkowitz & Johnson 1995) are identical to ours to lowest order in rj. Remarkably, though, their reported theoretical estimates 



of the system frequencies are not quite as accurate as ours (Merkowitz & Johnson 1997). The reason for this is probably the 



following: in M&J's model these frequencies appear within an algebraic system of 5+ J linear equations with as many unknowns 
which has to be solved; in our model the algebraic system has only J equations and unknowns, actually equations (^3|). This is 
a very appreciable difference for the range of values of J under consideration. While the roots for the frequencies can be seen 
to mathematically coincide in both approaches, in actual practice these roots are estimated, generally by means of computer 
programmes. It is here that problems most likely arise, for the numerical reliability of an algorithm to solve matrix equations 
normally decreases as the rank of the matrix increases. The significant algebraic simplification of our model's equations should 
therefore be considered one of real practical value. 

7.2 Other mismatched parameters 

We now assess the system sensitivity to small mismatches in resonators' masses, locations and frequencies. 
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Figure 6. The frequency spectrum of the TIGA distribution as resonators are progressively added from none to 6. Continuous lines 
correspond to measured values, and dashed lines correspond to their r;^/-^ theoretical estimates with equation (p2). 
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7.2.1 Resonator mass mismatches 

If the masses are slightly non-equal then we can write 

Ma = vM{l + fiaV^^^) , a = l,...,J (66) 

where rj can be defined e.g. as the ratio of the average resonator mass to the sphere's mass. It is immediately obvious from 
equation ( |6^ that mass non-uniformities of the resonators only affect our equations in second order, since resonator mass 
non-uniformities result, as we see, in corrections of order i]^^^ to r]^^^ itself, which is the very parameter of the perturbative 
expansions. The system is thus clearly robust to mismatches in the resonator masses of the type (pw. 



7.2.2 Errors in resonator locations 

The same happens if the locations of the resonators have tolerances relative to a preselected distribution. For let be a set 
of resonator locations, for example the TIGA or the PHC positions, and let n'a be the real ones, close to the former: 

n'a = ria -f Va 77^''^ , a = l,...,J (67) 



The values Ha determine the eigenvalues (^a in equation (|31|) , and also they appear as arguments to the spherical harmonics 
in the system response functions of sections ^|-|^. It follows from ( |6^ ) by continuity arguments that 

Yimin'a) = Yi,„{na) + 0{v^/^) (68) 
C = Ca + Oiv'^^) (69) 

Inspection of the equations of sections U-o shows that both ("a and Yim{na) always appear within lowest order terms, and 
hence that corrections to them of the type (|68f)-(p9[) will affect those terms in second order again. We thus conclude that the 
system is also robust to smaU misalignments of the resonators relative to pre-established positions. 



7.2.3 Resonator frequency mistunings 



The resonator frequencies may also differ amongst them, so let 

na = fi{l + PaV^^^) , a = l,...,J (70) 

To assess the consequences of this, however, we must go back to equation ([B6| ) and see what the coefficients in its series 
solutions of the type (BE) are. The procedure is very similar to that of section 7.1, and will not be repeated here; the lowest 
order coefficient Xi seen to satisfy the algebraic equation 



det 



Sab 



1 V "-i 

J_ Xa£. 



(71) 



which reduces to ( |B11[ ) when all the p's vanish, as expected. This appears to potentially have significant effects on our results 
to lowest order in 77^^, but a more careful consideration of the facts shows that it is probably unrealistic to think of such large 
tolerances in resonator manufacturing as implied by equation (^^ in the first place. In the TIGA experiment, for example 
(Merkowitz 1995), an error of order r^^"^ would amount to around 50 Hz of mistuning between resonators, an absurd figure 
by all means. In a full scale sphere (~40 tons, ~3 metres in diameter, -^800 Hz fundamental quadrupole frequency, 77 ~ 10~^) 
the same error would amount to between 5 Hz and 10 Hz in resonator mistunings for the lowest frequency. This is probably 
excessive for a capacitive transducer, but may be realistic for an inductive one. With this exception, it is thus more appropriate 
to consider that resonator mistunings are at least of order 77. If this is the case, though, we see once more that the system is 
quite insensitive to such mistunings. 

Summing up the results of this section, we can say that the resonator system dynamics is quite robust to small (of order 
77^''^) changes in its various parameters. The important exception is of course the effect of suspension drillings, which do result 
in significant changes relative to the ideally perfect device, but which can be relatively easily calculated. This theoretical 
picture is fully supported by experiment, as robustness in the parameters here considered has been reported in reference 
([Merkowitz fc Johnson 1997|). 



8 CONCLUSIONS 

A spherical GW antenna is a natural multimode device with very rich potential capabilities to detect GWs on earth. But 
such detector is not just a bare sphere, it requires a set of motion sensors to be practically useful. It appears that transducers 
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of the resonant type are the best suited ones for an efficient performance of the detector. Resonators however significantly 
interact with the sphere, and they affect in particular its frequency spectrum and vibration modes in a specific fashion, which 
has to be properly understood before reliable conclusions can be drawn from the system readout. 

The main objective of this paper has been the construction and development of an elaborate theoretical model to describe 
the joint dynamics of a solid elastic sphere and a set of radial motion resonators attached to its surface at arbitrary locations, 
with the purpose to make predictions of the system characteristics and response, in principle with arbitrary mathematical 
precision. 

We have shown that the solutions to our equations of motion are expressible as an ascending series in powers of the small 
"coupling constant" rj, the ratio of the average resonator mass to the mass of the large sphere. The lowest order app roximation 
corresponds to terms of order r/^^'^ and, to this order, we recover, and widely generalise, other authors' results (Merkowitz 
& Jolnson 1997; Stevenson 1997; Magalhaes et al. 1997), obtained by them on the basis of certain simplifying assumptions. 
This has in particular enabled us to assess the system response for arbitrary resonator layouts, and to search the equations for 
configurations other than the highly symmetric TIGA. This search has led us to make a specific proposal, the PHC, which is 
based on a pentagonally symmetric set of 5 rather than 6 resonators per quadruopole mode sensed. The PHC distribution has 
the very interesting property that mode channels can be constructed from the resonators' readouts, much in the same way as 
in the TIGA ( VIerkowitz fc Johnson 1995| ). In the PHC however a new and distinctive characteristic is present: different wave 
amplitudes selectively couple to different detector modes having different frequencies, so that the antenna's mode channels 
come at different rather than equal frequencies. The PHC philosophy can be extended to make a multifrequency system 
by using resonators tuned to the first two quadrupole harmonics of the sphere and to the first monopole, an altogether 11 
transducer set (Lobo & Serrano 199C). 

The assessment of symmetry failure effects, as well as other parameter departures form ideality, has also interested us 
here. This is seen to receive a particularly clear treatment in our general scheme: the theory transparently shows that the 
system is robust against relative disturbances of order r] or smaller in any system parameters, and provides a systematic 
procedure to assess larger tolerances — up to order The system is shown to still be robust to tolerances of this order in 

some of its parameters, whilst it is not to others. Included in the latter group is the effect of spherical symmetry breaking 
due to system suspension in the laboratory, which causes degeneracy lifting of the sphere's eigenfrequencies, now split up into 
multiplets. By using our algorithms we have succeeded in numerically reproducing the reportedly measured frequencies of 
the LSU prototype antenna ( Merkowitz 199f^ ) with a fully satisfactory precision of four decimal places. The experimentally 
reported robustness of the system to resonator mislocations (Merkowitz & Johnson 1997) is also in full agreement with our 
theoretical predictions. 

The perturbative approach we have adopted is naturally open to refined analysis of the system response in higher orders 
in ri. For example, we can systematically address the weaker coupling of non- quadrupole modes. It app ears however that such 
refinements will be largely masked by noise in a real system, as shown by Merkowitz and Johnson ( Merkowitz fc Johnson 
199^ and this m ust therefore be considered first. So our next step is to include noise in the model and see its effect. Stevenson 
( Stevenson 1997 ) has already made some progress in this direction, and partly assessed the characteristics of TIGA and PHC, 
but more needs to be done since not too high signal-to-noise ratios should realistically be be considered in an actual GW 
detector. In particular, the discovery of mode channels also for the PHC distribution opens the possibility of analysis of noise 
correlations and dependencies, as well as th e errors in GW parameter estimatio n. These are natural extensions of this research, 
and some of them are currently underway ( Merkowitz, Lobo fc Serrano 19991 ). 
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APPENDIX A: GREEN FUNCTIONS FOR THE MULTIPLE RESONATOR SYSTEM 

The density of forces in he rhs of equation (|^) happens to be of the separable type 

f(x,t) = Vf(°'(x)g(°'(t) (Al) 



where q is a suitable label. We recall from reference ( Lobo 1995 ) that, in such circumstances, a formal solution can be written 
down for equation (^) in terms of a Green function integral, whereby the following orthogonal series expansion obtains: 

"(^' = E E /ir ' (x) (t) ( A2) 

a N 

where 



1 



^ / u^(x).f^'^J(x)d^x (A3) 

•J Sphere 

g'^^Ht) = [ g^^\t')smcuMit^t')dt' (A4) 
Jo 

Here, ujn and U]v(x) are the eigenfrequencies and associated normalised wavefunctions of the free sphere. Also, is an 
abbreviation for a multiple index {nlm}. The generic index a is a label for the different pieces of interaction happening in 
the system. We quote the result of the calculations of the terms required in this paper: 



/rcsonators.iV = l^a- U% {^a)] , a=l,...,J (A5) 

/gw^jv = a„i5w5mm' , N = {nlm} , I' ^0,2, m' ^ -I' , . . . ,1' (A6) 

/stroke, JV = fo-ulfixo) (A7) 

© 1999 RAS, MNRAS 000, 



Multiple mode gravitational wave detection with a spherical antenna 21 

where the coefficients a„i in ( [A6[ ) are overlapping integrals of the type ( [A3| ) |* * *[ and 



o*"^ (t) 

resonators, A/" V / 



9gw,n 



9stTokc,N 



[za{t') — Ua(t)] sin UJ N{t — t') dt' 
g^^"^\t') sin LUN{t — t') dt' 



sin Art 



If this is replaced into one readily finds 
'' ' M, 



U(x,t) = "^Jv'uivW {Yj^^b [ni,-U;,(Xi,)] P^eLnators.ivW 



•/external, AT ^external, iV V^/ 



(A8) 

(A9) 
(AlO) 

(All) 



where the label "external" explicitly refers to agents acting upon the system from outside; two kinds of such external actions 
are considered in this article: those due to GWs and those due to a calibration hammer stroke. Specifying x = Xa in and 
multiply by iia above, the following is readily found: 



Uait) 



^external(^) + ^ r,, / K^tit - t') [z,{t') - U,(t')\ dt' 

b=i Jo 

where r,t = Mt/M, Mr"='""'(^)= n„-u°'''<=^"=''(x„ t), 

^externalj.^^ t) = J^Yl "^^^ /ixternal.iV "JV (x) !?external,]V (*) 
a N 

and 

Kab{t) = Ol ^ Lo'^^ [nb-u%{xt)] [na-U]v(xa)] siu CJjvt 



(A12) 
(A13) 

(AM) 
(A15) 



The following bare sphere responses to GWs and hammer strokes (equations (^ and (^o|), respectively) can be calculated 
by direct substitution. We present here the results as Laplace transform domain functions, which are the useful ones for our 
purposes: 



-GW/ N 

Ua (a) 



" stroke / \ 



E E 



1^0 and 2 \n—l 
m— — l,. . . ,/ 



Yimina) p*''"'(s) , a = 1, 



E 

nl 



fo 



\Ar,i{R)\^ Pi{na-no) , a = l,...,J 



(A16) 



(A17) 



where Yim are spherical harmonics and Pi Legendre polynomials ( Edmonds 1960| ). The calculation of the Laplace transform 
of the kernel matrix (A15) is likewise immediate: 



Kabis) = o^'' 9 [nt-ujf (Xi,)] [na-Ujv(Xa)] 



Given that (see ( Lobo 199E ) for full details) 
u„im (x) = A„i (r) Yijn{9,ip)n- B^i (r) in x LF;^ {9, ip) 



(A18) 



(A19) 



and that the spheroidal frequencies Lj„i are 2Z+l-fold degenerate, (A18) can be easily summed over the degeneracy index m, 
to obtain 



(A20) 



* * * We use here the better adapted notation a„o and a„2 instead of a„ and 6„, respectively, of reference (Lobo 1995). Numerical 
values are however identical. 
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or, equivalently, 



Katis) = J2 



S2 



47r 



Pi(na-n6) 



(A21) 



where use has been made of the summation formula for the spherical harmonics (Edmonds 196C) 



y;m(n6) Vi^(na) = — ^— P,(na-n6) 

m— — Z 

and where Pi is a Legendre polynomial: 



(A22) 



(A23) 



APPENDIX B: SYSTEM RESPONSE ALGEBRA 

From equation (^^, i.e., 



E 



+ v 



(ni) 
2 ^ ''^ab 



s2 + f72 



s2 + n2 ' 



- UJriglg) 



(Bl) 



we must first isolate qb{s), then find inverse Laplace transforms to revert to time domain quantities. Substituting the values 
of M?^(s) and ^(s) from ([aI6|) and Q into we find 



q4.)= ^ , a = l,...,J 



(B2) 



i = and 2 

771 — 



where 



S2 



-;2+E 



a77i A„i{R) 



E 



Sab + V ^ 



(s2 + J72)(s2 



_ ("0 

,2 -''^ab 



yi7n(ni,) 



(B3) 



Now, using the convolution theorem of Laplace transforms, we see that the time domain version of equation (B2) is 



?4i) = E f ^^a'"\t-t')g^'"^\t')dt' 
I. Jo 



a — 1, 



(B4) 



where <&i''"'(t) is the inverse Laplace transform of (|B>j). The inverse Laplace transform of $a (s) can be expediently 



(Im) I 



calculated by the residue theorem through the formula ( Helstrom 1968| ) 

= 2-Ki ^ ^ jresidues of (s) e"*] at its poles in complex s — plane} 



(B5) 



Clearly thus, the poles of ^J, '" (s) must be determined in the first place. It is immediately clear from equation (B3) that 
there are no poles at either s = 0, or s = ±iQ., or s = ±iu>ni, for there are exactly compensated infinities at these locations. 
The only possible poles lie at those values of s for which the matrix in square brackets in (B3) is not invertible, and these 
correspond to the zeroes of its determinant, i.e., 

s^ 

-^-]-^Kab{s) =0, poles (B6) 



A(s) = det 



5ab + 11 



There are infinitely many roots of equation (Bf]), and analytic expressions cannot be found for them. Perturbative 
approximations in terms of the small parameter ri will be applied instead. It is assumed that 

Q = Ulnolo (B7) 



for a fixed multipole harmonic {nolo}. Equation (B6) can then be recast in the more convenient form 



det 



Sab + V 



(s2 + f72): 



■ Xab 



+ n E 



(,2 + f^2)(,2+^2J 



(7l!) 



Xab 



= 



(B8) 



Since 77 is a small parameter, the denominators of the fractions in the different terms in square brackets in (BJ) must be 
quantities of order r) at the root locations for the determinant to vanish at them. A distinction however arises depending on 
whether s'^ is close to —Q? or to the other —ui^i- There are accordingly two categories of roots, more precisely: 
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si 



2 



1'' (l + XI v'^^ + Xiri+--) {^^ ^noio) (B9) 
li (l + b["'^ V + 4"'' + (nl / noio) (BIO) 

can be calculated recursively, starting form the first, by substitution of the 



The coefficients x^, Xiv a^nd fej"'' 
corresponding series expansions into equation (B5). The lowest order terms are easily seen to be given by 



det 

and 
det 



= 



—2 Ol Oab - Xab 



= 



(Bll) 



(B12) 



nl) 
b 



respectively. Both equations (Bff) and (Bf2) are algebraic eigenvalue equations. As shown in appendix |c|, the matrix xi' 

has at most {21 + 1) non-null positive eigenvalues — all the rest up to J are identically zero. 

As a final step we must evaluate (|B5| ). This is accomplished by standard textbook techniques (see e.g. ([Porter fc Stirling 
199C0; the algebra is quite straightforward but rather lengthy, and we shall not delve into its details here, but quote only the 



most interesting results, ft appears that the dominant contribution to ^a"^\t) comes from the poles at the locations ( [B9| ), 

(Bf3) 



whereas all other poles only contribute as higher order corrections; generically, $a"^\t) is seen to have the form 

0(0) 



'fa ""^ (i) oc 77 y ^ {sincuc+t — sinujc-t) S, 



llo 



where 



2 



I ± 



\Anoio{R)\ C^ri'A +o{n) , 



In Laplace domain one has, 

e"'(.)oc,7-^/^ ^ [ 



2 , 2 \ -1 
S +UJc+) 



/ 2 , 2 



(B14) 



(B15) 



Detailed calculation of the residues ( Serrano f QQE ) yield equation (^0|), which must be evaluated for each particular tuning 
and resonator distribution, as described in section ^ 



APPENDIX C: EIGENVALUE PROPERTIES 

We give in this Appendix a few important properties of the matrix Pi{na-nb) for arbitrary I and resonator locations ria 
(a=f ,. . . ,J) which are useful for detailed system resonance characterisation. 

We first note the summation formula for spherical harmonics (Edmonds f960) 

Y*^{iia)Yim{nb) = Pi{na-nb) , a, 6 = 1,..., J (CI) 

m — — l 

where Pi is a Legendre polynomial: 

^K^) = ^ ^ (^^ - D' (C2) 

To ease the notation we shall use the symbol Vi to mean the entire JxJ matrix Pi{na-ni,), and introduce Dirac kets \m) 
for the column J-vectors 

/ YUni) \ 

; , m^-l,...,l (C3) 

\ Yi,4nj) J 

These kets are not normalised; in terms of them equation ( |Cl|) can be rewritten in the more compact form 

m— — Z 
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Equation (C4) indicates that the rank of the matrix Vi cannot exceed {21 + 1), as there are only {21 + 1) kets |m). So, if 
J > {21 + 1) then it has at least {J — 21 — 1) identically null eigenvalues — there can be more if some of the iia's are parallel, 
as this causes rows (or columns) of Vi to be repeated. 

We now prove that the non-null eigenvalues of Vi are positive. Clearly, a regular eigenvector, |<^), say, of Vi will be a 
linear combination of the kets | m) : 

I 

m — — Z 

where we have called (^^ the corresponding eigenvalue, since it is a positive number, as we shall shortly prove. If the second 
( |C5|) is substituted into the first then it is immediately seen that 

I 

(C^'^mm' - ("ijm')) (^^/ = (C6) 

m' — — Z 

which admits non-trivial solutions if and only if 

det (C' (5,„™' - (m|m')) = (C7) 

In other words, are the eigenvalues of the {21 -I- 1) x {21 -\- 1) matrix (mjm'), which is positive definite because so is the 
"scalar product" {4>\4>). All of them are therefore strictly positive. 
Finally, since the trace is an invariant property of a matrix, and 
J ./ 

trace(Pi) P,(na-na) = ^ 1 = J (C8) 

a — 1 a — 1 

we see that the eigenvalues ("^ add up to J: 
./ 

trace(PO = = I] Ca = J (C9) 

a = l C„5^0 

This paper has been produced using the Royal Astronomical Society/Blackwell Science ffl^njX style file. 
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